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Hexagonal  Ferrite  Ku-Band  Radar  Absorbent 
Materials 

92FE0848A  Chengdu  DIANZI KEJI DAXUE  XUEBAO 
[JOURNAL  OF  UNIVERSITY  OF  ELECTRONIC 
SCIENCE  AND  TECHNOLOGY  (UEST)  OF  CHINA] 
in  Chinese  Vol  21  No  i,  Jun  92  pp  312-316 


[Article  by  Guo  Bijun  [6665  3880  0689],  Deng  Longjiang 
[6772  7893  3068],  and  Chen  Xuanzhen  [7115  6693 
5271]  of  the  Department  of  Materials  Science  and  Engi¬ 
neering,  UEST,  Chengdu  610054:  “Hexagonal  Ferrite 
Ku-Band  Radar  Absorbers”;  MS  received  8  Sep  91, 
revised  28  Oct  91;  project  supported  by  the  Preliminary 
Research  Fund  of  MMEI] 

[Text]  Abstract 

The  design  concept  and  preparation  technique  of  hexag- 
onal-ferrite-based  Ku-band  radar  absorbent  materials 
(RAM)  is  discussed.  The  physical  properties  and  mea¬ 
sured  electromagnetic  parameters  are  presented.  Finally, 
preliminary  application  results  are  described.  The  absor- 
bence  is  greater  than  12  dB  in  the  2-cm  band.  If  the 
composition  of  the  coating  is  further  adjusted,  an  absor- 
bence  of  greater  than  10  dB  can  be  achieved  in  the  2-3 
cm  band  (8-18  GHz)  range. 

Key  Words:  Ku-band,  ferrite,  absorber,  spontaneous 
resonance,  radar  absorber,  coprecipitation  technique 

Classification:  TM277 


I.  Introduction 

The  use  of  magnetic  materials  to  absorb  microwave  energy 
is  based  on  spontaneous  resonance.  When  the  microwave 
frequency  coincides  with  the  natural  resonance  frequency  of 
the  material,  maximum  absorption  occurs.  The  spontaneous 
resonance  frequency  is  proportional  to  the  anisotropic  mag¬ 
netic  field  inside  the  material.  Hexagonal  ferrites  have  a 
high  anisotropic  field.  Therefore,  their  natural  resonance 
frequencies  are  high  and  can  be  used  to  absorb  microwave 
energy.  Elements  such  as  Co  and  Ti  are  used  to  replace  the 
Fe^  ion  in  BaFeOi9  to  produce  a  microwave  absorber 
Ba(CoTi),tFei2.xOi9  substitution.  The  anisotropic 

field  can  vary  over  a  wide  range  as  a  function  of  x. 


Figure  1.  Ba(CoTi)j^Fei2.xOi9  Natural  Resonance  Fre¬ 
quency  vs.  X 


Consequently,  it  is  possible  to  produce  a  wide-band  micro- 
wave  absorber.  Figure  1  shows  how  the  natural  frequency  of 
Ba(CoTi)xFei2.xOi9  decreases  as  a  function  of  x. 

In  addition  to  composition,  the  shape,  size  and  orienta¬ 
tion  of  the  absorbent  power  have  significant  effects  on  its 
absorption  characteristics.  When  the  power  size  is 
reduced  to  the  monodomain  level,  the  domain-wall 
resonance  mechanism  no  longer  exists.  At  this  time,  the 
natural  resonance  frequency  of  the  power  increases 
slightly.  However,  the  maximum  microwave  energy  it  is 
capable  of  absorbing  falls  slightly  as  well,  as  shown  in 
Figure  2.  The  powder  prepared  in  this  work  is  approxi¬ 
mately  7-10  pm  in  size. 


Average  particle 
size 


Figure  2.  Effect  of  Powder  Size  on  Microwave  Absor- 
bence  a 


Orientation 
field  =  160  A/m 


0  IC  12  14  16  /(GHz) 


Figure  3.  Effect  of  Magnetic  Field  Orientation  on 
Microwave  Absorbence  a 


11.  Experimental  Results 

High-purity  BaC03,  Fe203,  CoO  and  Ti02  are  used  as 
raw  materials.  They  are  weighed  stoichiometrically  and 
then  ground  in  a  ball  mill  for  8  hours.  The  powder  is 
pre-heated  at  1280-1350°C,  crushed  and  then  sintered. 
After  it  is  crushed  again  and  sieved,  it  finally  goes 
through  an  annealing  process.  The  physical  properties 
and  electromagnetic  parameters  of  different  M-phase 
powders  of  different  x  content  were  experimentally 
determined.  The  results  are  as  follows: 
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1)  A  scanning  electron  micrograph  of  the  powder,  shown 
in  Figure  4,  clearly  proves  that  the  powder  is  hexagonal 
and  the  average  particle  size  is  a  <  10  pm. 


Figure  4.  Scanning  Electron  Micrograph  of  a 
Hexagonal  Ferrite  Ku-Band  Radar  Absorber 


2)  The  X-ray  diffraction  analysis  result,  shown  in  Figure  5, 
indicates  that  the  material  is  a  single  M-phase  material 

3)  The  specific  saturated  magnetization  moment  Og  is 
measured  using  a  magnetic  balance  and  Og  =  54.0  (emu/g). 

4)  An  HP85 1  OB  network  analyzer  was  used  to  determine 
electromagnetic  parameters  such  as  relative  complex 


Figure  5.  X-Ray  Diffraction  Pattern  of  Ferrite  Ku-Band 
Absorber 


magnetic  permeability  Pj.  =  p/  -  jPr“  and  relative  com¬ 
plex  dielectric  permeability  =  e/  The  results  are 
shown  in  Table  1. 


Table  1.  Relative  Complex  Magnetic  Permeability  and  Dielectric  Permeability,  HP8510B  Network  Analyzer 


F(Gllz) 

Er’ 

Er“ 

Ef^/Er’ 

w 

Ur“ 

Ur’VUr’ 

12.40 

6.955 

0.076 

0.011 

1.796 

0.323  i 

. 

0.100 

13.50 

7.032 

0.079 

0.011 

1.784 

0.727 

0.407 

14.50 

7.036 

0.121 

0.017 

1.360 

0.998 

0.729 

15.50 

7.031 

0.100 

0.015 

1.001 

0.830 

16.50 

7.021 

0.119 

0.017 

0.924 

0.616 

0.667 

17.50 

7.021 

0.119 

0.017 

0.954 

0.520 

0.546 

18.00 

7.078 

0.061 

0.007 

0.956 

0.516 

0.540 

III.  Application 

1)  The  initial  Ku-band  RAM  made  have  no  other  additives  besides  an  epoxy  binder.  Table  2  shows  the  reflective 
characteristics  of  the  coated  specimens. 


Table  2.  Reflective  Characteristics  of  RAM  Made  of  Ferrites  and  Epoxy  Binder 
(absorber  thickness  1.5  mm,  area  180  x  180  mm^  single  layer,  backed  by  metal  plate) 


Frequency  (GHz) 

Reflectivity  (dB) 

Frequency  (GHz) 

Reflectivity  (dB) 

12.40 

-3.22 

15.40 

-9.38 

12.60 

-3.86 

15.60 

-9.14 

12.80 

-4.62 

15.00 

-8.62 
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Table  2.  Reflective  Characteristics  of  RAM  Made  of  Ferrites  and  Epoxy  Binder 


(absorber  thickness  1.5  mm,  area  180  x  180  mm^  single  layer,  hacked  by  metal  plate)  (Continued) 


Frequency  (GHz) 

Reflectivity  (dB) 

Frequency  (GHz) 

Reflectivity  (dB) 

13.00 

-6.24 

16.00 

-8.06 

13.20 

-6.98 

16.20 

-8.48 

13.40 

-7.68 

16.40 

-8.74 

13.60 

-8.58 

16.60 

-8.90 

13.80 

-9.72 

16.80 

-9.34 

14.00 

-10.50 

17.00 

-9.52 

14.20 

-10.74 

17.20 

-9.90 

14.40 

-10.88 

17.40 

-10.06 

14.60 

-10.78 

17.60 

-9.86 

14.80 

-10.26 

17.80 

-10.54 

15.00 

-9.90 

18.00 

-4.04 

15.20 

-9.86 

2)  The  reflective  characteristics  of  specimens  coated  with  a  combination  of  ferrite  and  other  absorbers  are  shown  in 
Tables  3  and  4. 


Table  3.  Reflective  Characteristics  of  Microwave  RAM  Made  of  Ferrite  and  Other  Absorbers 
(absorber  thickness  1.29  mm,  area  180  x  180  mm^  single  layer,  backed  with  metal  plate) 


Frequency  (GHz) 

Reflectivity  (dB) 

Frequency  (GHz) 

Reflectivity  (dB) 

12.40 

-10.20 

15.40 

-20.42 

12.60 

-11.40 

15.60 

-22.52 

12.80 

-12.70 

15.80 

-22.26 

13.00 

-15.60 

16.00 

-24.92 

13.20 

-17.02 

16.20 

-24.54 

13.40 

-20.60 

16.40 

-24.66 

13.60 

-24.82 

16.60 

-24.42 

13.80 

-26.28 

16.80 

-22.88 

14.00 

-24.50 

17.00 

-21.82 

14.20 

-22.78 

17.20 

-19.52 

14.40 

-21.02 

17.40 

-17.64 

14.60 

-20.54 

17.60 

-16.10 

14.80 

-19.62 

17.80 

-14.50 

15.00 

-20.38 

18.00 

-10.78 

15.20 

-20.86 

Table  4.  Reflective  Characteristics  of  Microwave  RAM  Coated  With  Ferrites  and  Other  Thickening  Agents 
(absorber  thickness  1.26  mm,  area  180  x  180  mm^  single  layer,  hacked  with  metal  plate) _ 


Frequency  (GHz) 

Reflectivity  (dB) 

Frequency  (GHz) 

Reflectivity  (dB) 

12.40 

-12.06 

15.40 

-19.58 

12.60 

-12.98 

15.60 

-20.56 

12.80 

-13.80 

15.80 

-19.52 

13.00 

-15.68 

16.00 

-20.66 

13.20 

-16.94 

16.20 

-20.42 

13.40 

-17.90 

16.40 

-20.40 
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Table  4.  Reflective  Characteristics  of  Microwave  RAM  Coated  With  Ferrites  and  Other  Thickening  Agents 
(absorber  thickness  1.26  mm,  area  180  x  180  single  layer,  backed  with  metal  plate)  (Continued) 


Frequency  (GHz) 

Reflectivity  (dB) 

Frequency  (GHz) 

Reflectivity  (dB) 

13.60 

-19.12 

16.60 

-20.30 

13.80 

-19.72 

16.80 

-18.98 

14.00 

-19.20 

17.00 

-18.82 

14.20 

-19.34 

17.20 

-17.32 

14.40 

-18.56 

17.40 

■16.06 

14.60 

-18.94 

17.60 

-15.32 

14.80 

-18.82 

17.80 

-14.34 

15.00 

-19.58 

18.00 

-12.52 

15.20 

-20.14 

3)  Coated  specimens  weigh  approximately  3.52  kg/m^. 

From  the  reflective  characteristics  of  the  two  types  of 
coatings  shown  above,  we  can  see  that: 

1)  The  absorbence  of  the  specimen  can  reach  above  12 
dB  over  the  entire  2-cm  band.  This  indicates  that  such 
absorbers  can  cover  the  entire  Ku  band. 

2)  If  the  absorbent  coating  formula  and  electric  structure 
can  be  further  adjusted,  it  is  highly  probable  that  a 
wide-bandwidth,  high-absorbence  coating  which  covers 
the  2-3  cm  (8-1 8  GHz)  band  with  an  absorbence  of  more 
than  10  dB  can  be  made. 

3)  Compared  to  conventional  ferrite-based  RAM,  the 
Ku-band  coating  developed  is  thicker,  but  lighter  in 
weight.  This  is  especially  critical  for  an  aircraft. 

IV.  Conclusions 

Using  ion  substitution,  it  is  possible  to  control  the 
microwave  operating  frequency  by  replacing  some  of  the 
Fe^  in  BaFei20i9  The  new  hexagonal  ferrite  developed 
is  an  excellent  Ku-band  absorber.  The  successful  devel¬ 
opment  of  such  a  RAM  is  of  significant  importance  to 
the  preparation  of  wide-band  stealth  materials. 
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Analysis  of  Multilayer  Frequency-Selective 
Surfaces 

92FE0801F  Beijing  DIANZI XUEBAO 
[ACTA  ELECTRONICA  SINICA]  in  Chinese 
Vol  20  No  6,  Jun  92  pp  85-88 

[Article  by  Shen  Zhongxiang  [3088  1813  4382]  and  Hua 
Ron^i  [5478  2837  0823]  of  Nanjing  Institute  of  Aero¬ 
nautics:  “Analysis  of  Multilayer  Frequency-Selective 
Surfaces”;  MS  received  Nov  90,  revised  Aug  91] 

[Excerpts]  Abstract 

The  scattering  characteristics  of  multilayer  frequency- 
selective  surfaces  (FSS)  are  analyzed  by  the  method  of 
moments  in  conjunction  with  generalized  scattering 
matrix  theory.  Through  the  introduction  of  a  generalized 
waveguide,  a  unified  treatment  of  a  multilayer  FSS  with 
an  arbitrarily  shaped  element  is  realized.  As  examples, 
the  scattering  characteristics  of  several  FSS  structures 
are  computed.  The  results  are  in  good  agreement  with 
published  experimental  data. 

Key  Words:  electromagnetic  wave,  scattering,  frequency- 
selective  surface. 

I.  Introduction 

It  is  well  known  that  a  FSS  has  the  property  of  a 
bandpass  or  band-rejection  filter.  It  is  widely  used  in 
radomes,  reflector  antennas  and  beam-control  systems. 
A  great  deal  of  work  has  been  done  on  the  analysis  of 
single-layer  FSSs.  Many  methods  have  been  presented, 
including  mode  matching,  spectral  domain^  and  sub- 
domain  moment.^  However,  these  methods  are  universal 
and  have  their  limitations.  As  for  the  analysis  of  multi¬ 
layer  FSSs,  it  is  limited  to  two  layers.^’® 

In  this  work,  the  scattering  characteristics  of  a  multilayer 
FSS  are  analyzed  using  the  method  of  moments  and 
generalized  scattering  matrix  theory.  The  basic  idea  is  to 
introduce  the  generalized  waveguide  concept  and  con¬ 
sider  each  element  on  the  FSS  as  a  joint  for  two  gener¬ 
alized  waveguides.  The  intrinsic  mode  of  a  generalized 
waveguide  is  the  periodic  Floquet  mode.^^  The  method 
of  moments  is  used  to  find  the  generalized  scattering 
matrix  at  the  waveguide  joint.  Then,  the  scattering 
parameters  of  the  multilayer  FSS  can  be  calculated  based 
on  the  generalized  scattering  matrix. 
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II.  Theoretical  Analysis 

The  following  is  an  analysis  of  the  theory  using  an  infinite 
conductor  screen  with  a  periodic  two-dimensional  array  of 
holes.  For  a  two-dimensional  planar  bonded-plate  array, 
the  analysis  can  be  derived  from  the  results  of  this  work  on 
the  basis  of  duality.  The  structure  studied  is  shown  in 
Figure  1,  where  dj  and  d2  are  the  periods  along  the  x-axis 
and  s-axis,  respectively,  and  a  is  the  angle  between  the  two 
axes,  [passage  omitted] 


rectangular  arrays)  are  shown  in  Figure  5.  Other  results 
in  existing  literature  are  also  shown  in  the  figures.  It  is 
obvious  that  they  are  in  excellent  agreement. 


6  10  U  18  22 


//GHz 

Figure  3.  Transmission  Characteristics  of  Main  Mode 
on  Infinite  Plate  With  Two-Dimensional  Rectangular 
Array 

(di  =  7.6  mm,  d2  ^  15.2  mm,  a  =  90°,  a  =  1.27  mm,  b  » 
13.5  mm,  &yi  6^2  ==  1) 


III.  Numerical  Results 

The  theory  discussed  above  is  used  to  calculate  the 
reflection  and  transport  characteristics  of  a  plane  wave 
incident  upon  a  single-  and  a  double-layer  FSS.  Specifi¬ 
cally,  the  FSS  scattering  characteristics  of  infinite  plates 
with  an  array  of  rectangular  and  circular  holes  are 
calculated  in  this  work.  In  the  computation,  =  Mq. 
Furthermore,  in  order  to  be  universally  applicable,  the 
intrinsic  mode  of  the  cross  section  of  the  waveguide,  in 
the  shape  of  the  hole  on  the  plate,  is  chosen  to  express 
Mq.  In  the  actual  computation,  Q  =  10,  i.e.,  one  con¬ 
siders  the  lowest  10  modes  in  the  waveguides.  In  equa¬ 
tion  (1),  p  and  q  should  be  selected  so  that  tpq  <  1 1  k 
(where  k  is  the  wave  number  in  free  space).  Figures  3  and 
4  show  the  reflection  and  transmission  characteristics  of 
a  vertically  incident  TE  wave  on  a  single-layer  FSS.  The 
transmission  frequency  characteristics  of  a  vertically 
incident  TE  wave  on  a  double-layer  FSS  (comprised  of 


J _ 


t 


rf,/A 

Figure  5.  Power  Transmission  of  Double-Layer  Rectan¬ 
gular  Array  FSS 

(di  -  d2,  a  -  90°,  a  =  b  -  0.7  di,  t  =  0.5a,  =  £^2  “ 
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IV.  Conclusions 

The  scattering  characteristics  of  multilayer  frequency- 
selective  surfaces  (FSS)  are  analyzed  by  the  method  of 
moments  in  conjunction  with  generalized  scattering 
matrix  theory.  The  theory  can  be  applied  to  any  multi¬ 
layer  FSS  structure  with  any  shape  of  elements.  The 
analysis  done  in  this  work  is  based  on  the  assumption 
that  the  thickness  of  the  FSS  is  infinitely  thin.  Neverthe¬ 
less,  it  can  be  easily  extended  to  multilayer  FSSs  with  a 
finite  thickness. 
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[Text]  Abstract 

Motion  compensation  is  a  key  problem  in  ISAR  (inverse 
synthetic  aperture  radar).  Usually,  compensation  is 
made  using  an  isolated  point  on  the  target.  In  reality,  this 
isolated  point  may  not  exist.  This  paper  presents  a  new 
scheme  to  compensate  for  the  echo  of  a  continuous 
target. 

Key  Words:  inverse  synthetic  aperture  radar,  radar 
imaging,  motion  compensation, 

I.  Introduction 

ISAR  imaging  belongs  to  small-rotation-angle  observa¬ 
tion  imaging  with  a  rotation  angle  ranging  from  3°-5®.  It 
is  known*  that  targets  such  as  airplanes,  missiles  and 
ships  are  considered  as  a  series  of  discrete  scattering 
points  (perhaps  very  concentrated)  as  far  as  a  high- 
resolution  radar  is  concerned.  Within  an  azimuth  angle 
of  less  than  30°,  these  scattering  points  can  be  considered 
approximately  fixed  and  the  scattering  strength  is  a 
constant.  Hence,  ISAR  imaging  is  used  to  plot  the 
distribution  of  these  scattering  points.  If  the  target  is  on 
a  rotating  platform,  this  is  not  difficult  to  accomplish. 
Therefore,  that  case  will  not  be  discussed. 

A  real  moving  target,  as  far  as  a  radar  is  concerned,  has 
a  translational  component  in  addition  to  rotation  that 
contributes  to  imaging.  Hence,  it  is  necessary  to  select  a 
reference  point  on  the  target  prior  to  ISAR  imaging  in 
order  to  compensate  for  the  translational  component  of 
the  reference  point.  This  is  equivalent  to  placing  the 
reference  point  on  the  rotation  axis.  The  imaging  process 
is  carried  out  based  on  compensated  data. 

ISAR  imaging  can  handle  a  variety  of  signals.  In  this 
paper,  the  discussion  is  limited  to  narrow  pulses.  The 
results  can  be  extended  to  other  types  of  signals. 

II.  Concept  and  Purpose  of  Compensation 

In  ISAR  imaging,  there  is  relative  motion  between  target 
and  radar.  As  shown  in  Figure  1 ,  using  the  target  as  the 
reference,  i.e.,  choosing  a  point  on  the  target  as  a  still 
reference  point,  the  radar  moves  along  the  broken  line 
(with  dashes  and  dots)  shown  in  the  sketch.  As  far  as 
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rotating  platform  imaging  is  concerned,  the  radar  should 
be  in  circular  motion  along  the  arc  shown  in  the  sketch. 
Therefore,  the  purpose  of  motion  compensation  is  to 
provide  proper  adjustment  to  the  time  delay  of  the  echo. 
(Intensity  should  also  be  adjusted.  However,  if  relative 
distance  does  not  vary  much,  intensity  variation  can  be 
neglected.) 


Assume  an  echo  is  expressed  as  follows: 

s(t)  =  f(t)  exp(j(Oot)  (1) 

where  f(t)  is  the  complex  envelope  of  the  echo  and  coq  is 
the  angular  frequency  of  the  carrier. 

By  delaying  the  echo  by  a  time  interval  x,  the  delayed 
signal  is: 

s(t-t)  =  f(t-t)  exp(-jo)oT)  expCjWot)  (2) 


rotation).  However,  only  a  few  special  cases  have  iso¬ 
lated  scattering  points.  Many  targets  have  continuous 
echoes  that  are  a  superposition  of  echoes  from  numerous 
scattering  points. 

In  the  event  that  there  is  no  isolated  point,  a  point  from 
which  its  delay  and  phase  shift  can  be  determined  should 
be  chosen  as  the  reference  point.  In  the  following,  we 
shall  prove  that  it  is  feasible  to  use  the  “scattering  center 
of  gravity”  as  a  reference  point. 

III.  Calculation  of  Time  Delay  and  Phase  Shift  of 
Scattering  Center  of  Gravity 

The  longitudinal  and  transverse  coordinates  of  the  scat¬ 
tering  center  of  gravity  of  a  group  of  scattering  points  can 
be  calculated  based  on  their  positions  and  echo  intensi¬ 
ties.  In  reality,  only  echo  data  is  available.  The  problem 
is  how  to  calculate  the  coordinates  from  echo  data  alone. 

In  order  to  demonstrate  how  the  longitudinal  time  delay 
of  the  scattering  center  of  gravity  can  be  calculated  from 
echo  data,  let  us  use  the  case  shown  in  Figure  2  as  an 
example.  In  that  case,  the  echoes  do  not  overlap  on  the 
time  axis.  In  the  kth  observation,  if  the  echo  intensity 
and  time  delay  of  each  scattering  point  is  and  tj^  (i  = 
1,  2,  ...),  then  the  longitudinal  coordinate  T,,  of  its 
scattering  center  of  gravity  is 

(3) 


Usually,  radar  signals  are  processed  at  the  fundamental 
frequency  after  mixing.  After  multiplying  the  above 
equation  by  exp(-ja)ot),  we  obtain  the  delayed  complex 
envelope  f(t-T)  exp(-jcDo'c).  When  the  time  delay  x  is 
changed,  the  complex  envelope  f(t-x)  and  initial  phase 
factor  exp(-jcoot)  will  also  change.  However,  the  rate  of 
change  is  different.  For  an  X-band  radar,  a  delay  of  the 
order  of  a  millimeter  causes  significant  change  in  its 
initial  phase.  On  the  other  hand,  it  almost  has  no  effect 
on  the  complex  envelope  when  the  pulse  width  is  several 
dozen  centimeters.  Based  on  this  reason,  and  the  fact 
that  initial  phase  changes  come  around  every  In,  motion 
compensation  can  be  performed  in  two  steps.  The  first 
step  is  to  carry  out  time  delay  translational  correction 
(rough  compensation).  The  accuracy  is  of  the  order  of  a 
centimeter.  The  next  step  is  to  further  adjust  the  initial 
phase  of  the  corrected  echo  (fine  compensation).  The 
error  should  be  far  less  than  k/2.  Upon  completion  of 
these  two  steps,  the  initial  phase  is  essentially  accurate. 
However,  the  complex  envelope  delay  might  have  an 
error  of  the  order  of  a  centimeter.  This  should  not  have 
any  significant  effect  on  imaging. 

Hence,  if  the  target  has  an  echo  from  an  apparent 
isolated  scattering  point,  motion  compensation  is  not 
difficult.  It  can  be  done  by  simply  choosing  that  isolated 
point  as  a  reference  point  to  align  the  envelope  and  to 
change  the  initial  phase  so  that  the  phased  difference 
across  a  cycle  is  zero  (which  is  equivalent  to  the  axis  of 


Electromagnetic  waves 


Figure  2.  Target  Scattering  Points  and  Echoes 


However,  it  is  not  so  simple  in  real  life.  The  echoes  from 
the  scattering  points  of  a  continuous  target  do  overlap.  In 
this  case,  the  kth  echo  can  be  expressed  as: 

(4) 


where  v(t)  is  the  transmitted  signal.  The  complex  enve¬ 
lope  of  the  echo  is 
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where  a(t)  is  the  complex  envelope  of  v(t).  It  is  obvious 
that 

* 


•a*(t-t,*)exp(-/0..n)  (6) 


The  first  term  in  equation  (6)  is  the  power  of  each 
scattering  point.  The  second  term  represents  the  cross 
term.  In  practice,  it  is  only  necessary  to  consider  the 
overlapping  part  of  the  echoes.  As  mentioned  earlier,  tji^ 
varies  only  slightly  for  neighboring  echoes  and  has  no 
effect  on  a(t).  However,  phase  may  vary  considerably  in 
a  random  fashion.  In  an  approximation,  ®ilk  can  be 
treated  as  a  symmetrically  and  randomly  distributed 
variable  with  a  zero  mean  value  for  each  observation. 
Hence,  when  taking  the  statistical  average  of  (6),  the 
cross  term  is  zero,  i.e., 

£:i/*a)|^3=:aa,|a(t-t,*)|*  (7) 

f 


This  suggests  that  there  is  no  difference  in  the  calculation 
of  echo  power  compared  to  the  case  where  all  the  echoes 
do  not  overlap  from  the  standpoint  of  statistical  average. 
Nevertheless,  in  order  to  calculate  the  scattering  center 
of  gravity  (as  shown  in  equation  3),  it  is  necessary  to  use 

^  jo’.- 1 s k. 


The  echoes  from  the  target  can  be  divided  into  range 
rings  according  to  different  side  lobes.  There  might  be 
many  scattering  points  inside  a  range  ring  (with  a  spe* 
cific  til,  value).  During  continuous  observation,  due  to 
phase  variation,  the  echo  also  fluctuates.  However,  as  a 
statistical  average,  it  is  identical  to  the  situation  in  which 
every  range  ring  has  a  single  scattering  point.  In  other 
words,  when  calculating  the  time  delay  of  the  longitu¬ 
dinal  coordinate  of  the  scattering  center  of  gravity  with 
equation  (3)  based  on  the  echoes  obtained,  it  also  fluc¬ 
tuates.  However,  its  average  is  the  correct  value. 

As  a  verification,  we  designed  a  moving  target  with  20 
scattering  points.  Figure  3  shows  the  time  delay  of  the 
scattering  center  of  gravity  calculated  from  simulated 
echoes.  The  disturbance  is  apparent.  After  smoothing,  it 
is  very  close  to  the  original  design  curve.  We  also  added 


some  noise  to  the  simulated  echoes.  The  effect  of  noise 
amplifies  the  disturbance,  which  can  be  reduced  by 
smoothing. 


Figure  3.  Translational  Trajectory  of  Scattering  Center 
of  Gravity 


After  delay  compensation  of  the  complex  envelope,  the 
next  step  is  initial  phase  compensation.  The  motion  of  a 
target  can  be  divided  into  translational  and  rotational. 
The  periodic  phase  shift  of  the  echo  from  each  scattering 
point  is  the  sum  of  the  phase  shifts  produced  by  these 
two  types  of  motion.  The  translational  phase  shift  is  the 
same  for  every  scattering  point,  while  the  rotational 
phase  shift  mi^t  be  different  from  point  to  point.  With 
respect  to  the  rotation  axis  (reference  point),  rotational 
phase  shift  is  zero.  Therefore,  when  the  scattering  center 
of  gravity  is  chosen  as  the  reference  point,  the  phase  shift 
at  the  center  of  gravity  (not  necessarily  a  scattering 
point)  should  be  calculated  and  used  for  initial  phase 
compensation. 

If  rotation  alone  is  taken  into  consideration  and  a 
stationary  reference  point  is  chosen  as  the  origin,  then 
the  longitudinal  displacement  of  scattering  point  i,  Ay, 
generated  as  a  result  of  rotation  by  an  angle  Aep  is  (see 
Figure  4): 

Ayj  -  Ri  cos  ((pi  -  Aep)  -  Rj  cos  <Pi  =  yi  (cos  A(p  -  1)  +  Xj  sin 
A(p  (8) 


Figure  4.  Longitudinal  Displacement  of  Scattering  Point 
Caused  by  Target  Rotation 
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In  reality,  the  target  is  several  hundred  kilometers  away 
and  two  pulses  are  separated  at  most  by  a  few  millisec¬ 
onds;  the  A(p  generated  by  target  movement  within  this 
time  interval  is  usually  less  than  0.01°.  If  Xj  and  yj  are  10 
meters,  then  the  first  and  last  terms  of  equation  (8)  are  of 
the  order  of  a  micron  and  millimeter,  respectively. 
Hence, 

Ayj  »  X;  A<p  (9) 

As  far  as  the  ith  echo  is  concerned,  the  phase-angle 
difference  between  two  consecutive  observations  due  to 
target  motion  is 

d>i  =  (4jt/Xo)  Xi  A<p  (1 0) 

If  the  transverse  coordinate  of  each  scattering  point  is  Xy 
(i  =  1,  2...),  then  the  transverse  scattering  center  of 
gravity  is 


x,=  '^a,x,/^a,  (11) 


By  multiplying  (47i/Xo)  A(p  to  both  sides  of  the  equation 
above,  and  by  taking  equation  (10)  into  consideration,  it 
is  possible  to  obtain  the  phase-angle  difference  of  the 
transverse  scattering  center  of  gravity  after  the  target 
rotates  by  an  angle  of  Acp. 

(12) 


In  reality,  the  distribution  of  scattering  points  is 
unknown.  We  need  to  find  the  waveform 

observed.  If  it  is  simple  as  is  the  case  shown  in  Figure  2 
where  the  echoes  do  not  overlap,  the  phase  shift  Oi  can 
be  calculated  from  the  phase  change  of  consecutive 
echoes.  Then,  Oq  can  be  obtained  from  equation  (12) 
based  on  the  intensity  of  the  echo  Oj. 

If  the  scattering  center  of  gravity  is  chosen  as  a  stationary 
reference  point,  then  the  phase-angle  difference  between 
observations,  4>o,  should  be  zero.  Using  the  Oq  obtained 
from  equation  (12)  as  the  initial  phased  compensation 
should  produce  the  desired  result. 

In  real  life,  many  sub-echoes  may  overlap  at  the  same 
time.  In  this  case,  the  phase  changes  of  these  sub-echoes 
cannot  be  determined.  In  order  to  obtain  Oq  from  (12), 
the  easiest  way  is  to  replace  Oj  with  echo  power  at  a 
specific  time  (usually  in  an  interval  of  a  pulse  width)  and 
to  use  the  phase  difference  of  the  echo  at  this  time 
interval  between  two  consecutive  observations  as  3>i. 
Does  this  arrangement  make  sense?  Let  us  take  a  look  at 
the  case  where  only  the  echoes  from  two  scattering  points 
overlap. 

As  we  mentioned  several  times  earlier,  in  a  few  neigh¬ 
boring  observation  periods,  as  an  approximation,  the 
complex  envelope  can  be  considered  as  stationary.  This 


means  that  the  echo  vector  from  any  scattering  point  at 
any  instance  remains  unchanged.  Only  the  phase  angle 
changes.  Figure  5  shows  the  situation  involving  two 
scattering  points  in  two  consecutive  observations.  In  the 
figure,  ri  and  T2  are  the  amplitudes  of  the  echoes  at  those 
two  scattering  points  at  that  instance;  Rj  and  R2  pe  the 
synthetic  amplitudes  of  the  two  observations.  Using  the 
vector  of  ri  as  a  standard,  the  phase  angle  varies  by  A\|/. 
The  vector  angles  of  the  two  sub-echoes  in  the  two 
observations  are  0  and  0  +  A0,  with  0  =  02  -  0i  and  A0  = 
A02  -  A01,  where  0i,  02  and  A0i,  A02  are  the  phase  angles 
of  the  sub-echoes  and  the  phase-angle  differences  in  the 
second  observation,  respectively.  (The  variation  of  0i  is 
not  shown  in  the  figure.)  If  A0  and  Ay  are  both  small, 
then  geometrically  we  can  determine  that 

Ri  R2  «  R  (13) 

A\|;  ==  r2A0  cos  (0  -  x|f)/R  (14) 

If  Oj  and  Oj  in  equation  (12)  are  directly  replaced  by  the 
echo,  the  contribution  to  equation  (1 2)  at  that  moment  is 
R^(A\};  +  A0i).  After  substituting  it  into  equation  (14)  and 
taking  the  geometry  shown  in  Figure  5  into  account,  we 
have 

R^(Acp  +  A0i)  =  (r2^  +  rir2  cos  0)  A0  +  (r2^  +  ri^  +  2rir2  cos 
0)A0i  (15) 

Taking  the  average  of  the  above  equation  over  0,  we  get 
E[R2(A<p  +  A0i)]  =  ri^  A01  +  A02  (16) 

This  demonstrates  that  the  phase-angle  difference  of  the 


Figure  5.  Variation  of  Echo  From  a  Rotating  Target 


transverse  scattering  center  of  gravity  derived  from  the 
power  and  phase  difference  of  the  measured  waveform  is 
consistent,  from  the  standpoint  of  statistical  average, 
with  that  calculated  based  on  all  scattering  points.  There¬ 
fore,  the  phase  compensation  calculated  from  the  wave¬ 
form  must  be  filtered  to  remove  disturbance  from  it. 

The  above  is  the  situation  when  two  sub-echoes  overlap. 
It  is  easy  to  extend  to  more  overlapping  sub-echoes.  In 
the  above  discussion,  it  is  also  assumed  that  A0  ^  0. 
When  A0  has  a  finite  value,  direct  analysis  is  more 
difficult.  Through  computer  simulation,  it  was  found 
that  there  is  a  difference  between  the  statistical  average 
of  the  calculated  transverse  scattering  center  of  gravity 
and  the  real  value.  However,  as  long  as  the  absolute 
value  of  A0  <  0.5  rad  (which  covers  most  cases),  the 
error  is  less  than  5  percent.  This  kind  of  error  is  tolerable 
in  ISAR  imaging. 
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With  a  repetition  frequency  of  the  order  of  a  few 
milliseconds,  the  phase-angle  differences  of  various  scat¬ 
tering  points  caused  by  rotation  are  small.  However, 
translational  motion  will  cause  a  larger  but  equal  phase 
shift  at  every  scattering  point.  Hence,  the  phase-angle 
differences  of  a  scattering  point  within  a  period  are 
distributed  over  a  small  angular  range  and  this  angular 
range  is  randomly  distributed  over  [0,  In],  Since  phase 
angle  is  a  cyclic  function  of  2k,  the  angular  range  cannot 
have  any  2n  jumps  when  we  calculate  the  phase  differ¬ 
ence  of  the  scattering  center  of  gravity. 

rV.  Experimental  Results 

In  order  to  confirm  the  above  analysis,  a  great  deal  of 
simulation  and  experimental  work  was  done.  In  com¬ 
puter  simulation,  a  flying  scattering  point  model  and  an 
electromagnetic  scattering  model  were  employed.  Exper¬ 
imentally,  measured  data  from  a  rotating  model  was 
used.  Furthermore,  simulated  flight  data  was  also 
studied.  The  results  are  excellent.  Nevertheless,  we  are 
most  concerned  about  imaging  of  flying  targets. 

We  were  fortunate  to  obtain  a  set  of  Boeing  727  flight 
data  from  Professor  Steinberg  of  the  University  of  Penn¬ 
sylvania.  Figure  6  shows  the  image  generated  after 
motion  compensation  (aircraft  and  radar  parameters  are 
shown  in  the  figure).  Comparing  to  the  actual  silhouette 
of  the  aircraft  at  the  bottom,  it  is  not  difficult  to  come  to 
the  conclusion  that  the  imaging  effect  is  satisfactory. 
This  illustrates  that  motion  compensation  based  on  the 
scattering  center  of  gravity  is  a  viable  approach.  This 
technique  will  be  verified  with  more  experimental  data, 
especially  actual  moving  target  data,  in  the  near  future. 


(a)  Sketch  of  the  aircraft 


(b)  Aircraft  image 

Figure  6.  Imaging  Based  on  Motion  Compensated  Data 


It  should  be  pointed  out  that  the  center  of  mass  of  the 
echo  and  the  actual  scattering  center  of  gravity  coincide 


only  in  a  statistical  average  sense.  Hence,  it  is  necessary 
to  smooth  out  the  results  while  performing  complex 
envelope  delay  compensation.  This  technique  has 
proven  to  be  successful  in  a  variety  of  experimental 
conditions.  As  for  smoothing  of  initial  phase  compensa¬ 
tion,  computer  modeling  was  also  successful.  However, 
it  may  not  work  on  actual  data  because  there  are  other 
factors  affecting  the  initial  phase.  For  instance,  insta¬ 
bility  of  emission  frequency  will  cause  random  distur¬ 
bance  of  initial  phase  and  this  should  not  be  smoothed 
out. 

Motion  compensation  in  ISAR  imaging  needs  further 
research.  The  method  of  using  isolated  points  as  refer¬ 
ence  points  is  effective  when  there  are  isolated  scattering 
points.  However,  it  becomes  useless  with  targets  that 
have  densely  distributed  scattering  points.  In  this  case, 
the  scattering  center  of  gravity  method  has  an  edge. 
Furthermore,  the  more  scattering  points  there  are,  the 
more  even  the  distribution  becomes  and  the  more  favor¬ 
able  the  compensation  is.  In  other  words,  isolated  point 
and  scattering  center  of  gravity  complement  each  other. 
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[Excerpts]  Abstract 

When  there  is  a  relatively  isolated  and  strong  scatterer  in 
the  target,  ISAR  [inverse  synthetic  aperture  radar] 
motion  compensation  can  be  transformed  into  the  align¬ 
ment  of  each  return  and  the  phase  compensation  of  each 
return.  In  this  work,  the  return  phase  is  used  as  an 
approximation  to  estimate  the  return  Doppler  fre¬ 
quency.  In  addition,  a  novel  method  is  presented  to 
eliminate  any  ambiguities  of  Doppler  frequency.  The 
target  trajectory  is  re-established  using  the  unambiguous 
Doppler  frequency  to  align  the  range  of  each  return. 

Key  Words:  ISAR  ima^ng,  motion  compensation,  Dop¬ 
pler  frequency  ambiguity  resolution. 

L  Introduction 

In  microwave  imaging  with  a  phased  array  antenna,  the 
antenna  must  be  aimed  at  the  target  or  target  area  for 
focusing.  This  is  to  perform  phase  compensation  for  each 
element  and  then  find  the  coherent  sum.  SAR  (synthetic 
aperture  radar)  or  ISAR  imaging  can  be  treated  as 


imaging  with  a  large-aperture  phased  array  antenna. 
Therefore,  each  return  should  be  compensated.  The 
difference  is  that  because  the  aperture  is  extremely  large, 
it  is  not  only  necessary  to  perform  phase  compensation 
but  also  range  alignment  for  each  return.  The  combina¬ 
tion  is  considered  as  motion  compensation  in  this  paper. 

In  studying  ISAR  imaging,  because  it  is  difficult  to 
accurately  measure  the  trajectory  of  the  target,  motion 
compensation  has  always  been  a  challenge.  In  order  to 
solve  this  problem,  many  assumptions  have  been  postu¬ 
lated.  For  instance,  the  target  has  an  isolated  and  strong 
scatterer  and  this  small  scatterer  can  be  resolved  from 
the  longitudinal  range  of  the  wide  transmission  signal. 
Furthermore,  the  scatterer  has  a  small  transverse  dimen¬ 
sion  (which  is  less  than  the  required  transverse  resolu¬ 
tion).  Then,  the  small  scatterer  can  be  treated  as  a 
relatively  isolated  source  (a  point  scatterer  in  an  ideal 
situation).  It  has  been  demonstrated  in  reference  1  by 
means  of  statistical  analysis  of  a  great  deal  of  experi¬ 
mental  data  that  the  probability  that  such  a  relatively 
isolated  scatterer  exists  is  80  percent  when  the  range 
resolution  is  5r  =  1  m. 

For  targets  such  as  airplanes,  the  isolated  scatterers 
might  be  wing  tips,  engines  and  tails.  Under  the  assump¬ 
tion  that  such  an  isolated  scatterer  could  be  found,  many 
ISAR  imaging  experiments  on  real  flying  targets  or 
airplane  models  have  been  reported. 

Motion  compensation  in  ISAR  imaging  primarily 
involves  range  alignment  for  every  return.  A  spatial 
cross-correlation  method  and  frequency-domain  method 
were  used  for  range  alignment  in  reference  2.  Both 
methods  involve  a  great  deal  of  computation.  Here,  we 
will  introduce  a  simple  and  direct  method  for  range 
alignment  by  estimating  the  unambiguous  Doppler  fre¬ 
quency  in  order  to  re-establish  the  trajectory  of  the 
moving  target. 

Since  the  radial  velocity  of  the  target  is  not  constant,  i.e., 
its  Doppler  frequency  is  not  invariant,  the  Doppler 
frequency  estimated  based  on  the  returns  accumulated 
over  a  period  At  in  the  past  cannot  accurately  reflect  the 
Doppler  frequency  at  the  present  time.  In  this  work,  we 
use  the  phase  of  the  target  return  to  estimate  the  Doppler 
frequency  of  the  echo  and  present  a  new  method  to 
eliminate  Doppler  ambiguities. 

The  range  alignment  method  is  described  below. 

11.  Motion  Compensation 

In  the  following  discussion,  let  us  assume  that  the 
operating  wavelength  of  the  radar  isX  =  0.05  m,  and  the 
range  resolution  is  8r  =  1  m  for  the  1 50  MHz  wide-band 
signal  emitted  by  the  radar.  After  taking  the  computation 
load  and  the  transverse  resolution  required,  the  pulse 
repetition  frequency  (prf)  chosen  is  F^  =  100  Hz. 

By  using  a  coherent  measurement  radar,  we  can  obtain 
data  such  as  the  target  range  R’(t)  and  complex  target 
return  a’(t).  If  discrete  sampling  is  employed,  we  get: 
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Target  range  R’(i),  complex  target  return  a’(ij)  (where  i  = 

1,  2 . N  represents  the  ith  transmission  (or  reception), 

j  =  1 ,  2, M  represents  the  jth  range  sampling). 

a’(ij)  can  be  used  to  determine  the  mode  phase  v|f’(ij)  of 
the  target  return  phase  t|/(i  j),  where  (-2n  <  \|/’  <  0). 

When  there  is  a  relatively  isolated  scatterer  in  the  target 
(i.e.,  an  isolated  point),  it  is  assumed  that  the  target  range 
measured  corresponds  to  the  range  of  the  isolated  point. 
Assuming  that  random  error  in  range  measurement  is 
totally  due  to  thermal  noise,  then  the  error  eR(i)  obeys  a 
normal  distribution  N(0,Or),  where  Or  is  the  rms  (root 
mean  square)  value  of  the  normal  distribution.  When 
2or  =  5r/2,  i.e..  Or  =  6r/4,  the  probability  of  aligning  the 
return  to  the  same  range  unit  5r  is  95.45  percent.® 
Because  the  accuracy  of  ranging  cannot  meet  the  above 
requirement,  it  is  necessary  to  perform  range  alignment 
for  each  return,  [passage  omitted] 

III.  Computer  Simulation  and  Analysis  of  Results 

With  reference  to  a  six-point  model,  as  shown  in  Figure 
1,  computer  simulation  was  done  using  the  method 
described  above.  The  ratio  of  points  (1)  through  (6)  is 
1:1:1:0.5:0.6:0.8.  Each  point  is  rigidly  connected  and  the 
points  are  equally  spaced.  Let  us  assume  that  the  target 
height  is  Yq  “  8,000  m  and  the  vertical  distance  is  Zq  = 
-39.192  km,  then,  the  nearest  range  is  Rq  =  40  km,  as 
shown  in  Figure  2.  Assuming  that  the  target  is  moving 
linearly  from  point  A  to  B  at  a  constant  speed  of  V,  =  250 
m/s,  as  shown  in  Figure  3,  then  AB  =  ViT  =  250  m/s  x  4 
s  =  1,000  m.  The  transverse  resolution  is 

2^9  '^2AB/Ro'^^  (22) 


Here,  the  signal-to-noise  ratio  (SNR)  is  defined  as  the 
average  SNR  of  the  element  of  return  with  the  strongest 
average  power: 

SNRA=max|^;|]|fl(i,/)lVc^S}  (23) 


where  1  :<  j  :<  M  =  64,  N  =  400  and  is  the  standard 
deviation  of  noise. 

Figures  4-6  show  the  results  obtained  at  different  SNRs. 
The  result  shown  in  Figure  6  is  considerably  worse  than 
those  shown  in  Figures  4  and  5.  Furthermore,  there  is 
target  movement  in  the  longitudinal  direction.  This  is 
because  of  an  error  in  the  expansion  of  the  ambiguous 
Doppler  frequency  at  lower  SNRs  which  causes  the 
linearly  varying  Doppler  frequency  to  have  an  abrupt 
discontinuity  (up  or  down),  as  shown  in  Figure  7.  The 
dotted  line  in  the  figure  represents  the  original  linear 
Doppler  frequency  and  the  solid  line  indicates  the  dis¬ 
continuity.  This  causes  a  relatively  large  error  when  we 
reconstruct  the  range  and  consequently  produces  more 


Figure  1.  Target  Model 


Figure  2.  Relative  Position  Between 
Radar  and  Target 


Figure  3.  Target  Trajectory 


error  in  range  alignment.  Of  course,  this  makes  it  impos¬ 
sible  to  produce  an  accurate  image  of  the  moving  target. 
As  for  the  cause  of  the  error  in  Doppler  frequency 
expansion,  it  is  due  to  the  fact  that  the  noise-to-signal 
ratio  is  too  high  in  some  instances.  There  is  a  large  error 
in  the  ambiguous  Doppler  frequency  derived  from  the 
phase  of  the  return  which  leads  to  an  erroneous  Doppler 
frequency  expansion.  Further  study  on  this  issue  and 
more  effective  expansion  of  Doppler  frequency  will  be 
discussed  in  a  later  article. 

On  the  basis  of  the  above  discussion,  in  order  jo  r^^pn- 
struct  the  trajectory  of  a  moving  target,  the  unambiguous 
Doppler  frequency  should  be  accurately  estimated,  first. 
When  the  target  contains  a  relatively  isolated  and  strong 
scatterer,  the  key  issues  in  ISAR  imaging  are  how  to 
perform  range  alignment  for  every  return  wave  and  how 
to  find  the  range  element  for  the  isolated  scatterer.  , 
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IVe  Conclusions 

Under  the  premise  that  the  target  has  a  relatively  iso¬ 
lated  and  intense  scatterer,  this  paper  presents  a  way  to 
estimate  the  Doppler  frequency  of  the  return  by  its  phase 
and  introduces  a  novel  method  to  resolve  any  Doppler 
ambiguities.  Unambiguous  Doppler  frequency  is  then 
used  to  reconstruct  the  trajectory  of  the  moving  target  to 
perform  range  alignment  for  every  return.  Furthermore, 
the  isolated  scatterer  is  located  based  on  the  least  square 


of  normalized  return  amplitudes.  The  phase  of  the 
isolated  scatterer  is  then  used  as  the  standard  to  com¬ 
pensate  for  the  phase  of  the  return  from  each  range 
element.  This  method  has  been  found  to  be  effective 
through  theoretical  analysis  and  computer  simulation. 
However,  when  such  an  isolated  scatterer  does  not  exist, 
or  even  when  such  an  isolated  point  exists  but  the 
required  transverse  resolution  is  high,  this  phase  com¬ 
pensation  method  is  no  longer  effective.  This  problem 
will  be  discussed  in  future  articles. 
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[Text]  Abstract 

It  is  well  known  that  the  range  resolution  and  cross-range 
resolution  of  a  range-Doppler  imaging  radar  are  deter¬ 
mined  by  the  effective  bandwidth  of  the  transmitted 
signal  and  the  angle  of  rotation  of  the  object  relative  to 
the  line-of-sight  of  the  radar  during  the  coherent  pro¬ 
cessing  period.  The  purpose  of  this  work  is  to  lift  the 
limitation  of  conventional  FFT  [fast  Fourier  transform] 
range-Doppler  processing  to  significantly  improve  the 
resolution  of  a  range-Doppler  imaging  radar.  The 
approach  adopted  is  to  use  linear  prediction  for  data 
extrapolation  before  performing  Fourier  reconstruction. 
Preliminary  results  of  rotating  platform  imaging 
obtained  in  a  microwave  anechoic  chamber  indicate  that 
this  is  a  promising  technique  to  enhance  the  resolution  of 
range-Doppler  imaging  radar. 

Key  Words:  radar  imaging,  superresolution  processing, 
linear  prediction,  data  extrapolation. 


I.  Introduction 

There  are  a  variety  of  range-Doppler  radars  for  different 
applications.  Synthetic  aperture  radar  (SAR)  mounted 
on  a  moving  platform  is  used  in  air-to-ground  imaging. 
Inverse  synthetic  aperture  radar  (ISAR)  fixed  on  the 
ground  is  used  in  the  imaging  of  airplanes,  space  objects 
and  heavenly  bodies.  A  rotating  platform  radar  is  also 
fixed  on  the  ground.  It  can  also  be  used  in  imaging 
targets  on  the  rotating  platform.  These  range-Doppler 
radars  have  a  common  critical  parameter,  i.e.,  resolution 
of  image  reconstruction. 

A  range-Doppler  radar  gets  the  required  resolution  from 
the  range-Doppler  principle.  ‘  The  range  resolution  of  the 
image,  p^,  is  dependent  upon  the  effective  bandwidth  of 
the  transmitted  signal  B: 

p,  =  C/2B  (1) 

where  C  is  the  speed  of  li^t.  On  the  other  hand,  the 
cross-range  resolution,  p^,  is  dependent  upon  the  angle 
A0  through  which  the  object  rotates  relative  to  the  radar 
line-of-si^t  (RLOS)  in  the  coherent  processing  period. 
The  relation  between  p^.  and  A0  is  as  follows: 

Pe  =  X/2A0  (2) 


where  X  is  the  radar  wavelength.  It  should  be  pointed  out 
that  equations  ( 1 )  and  (2)  are  rigorously  valid  only  when 
the  A0  is  so  small  that  the  scattering  point  of  the  target 
does  not  go  over  a  range  resolution  unit  and  a  Doppler 
resolution  unit.  With  classic  signal  processing  tech¬ 
niques,  it  is  impossible  to  overcome  this  limitation 
imposed  by  equations  (1)  and  (2).  In  this  case,  the 
effective  bandwidth  of  the  transmitted  signal  and  the 
total  angle  of  rotation  of  the  target  with  respect  to  the 
RLOS  must  be  increased  to  enhance  the  range  resolution 
and  cross-range  resolution  of  the  image,  respectively.  It 
is  often  not  desirable  to  do  so  in  many  applications. 
Another  approach  to  improving  the  resolution  capability 
is  to  use  advanced,  rather  than  conventional,  signal 
processing  methods  because  these  methods  may  be  able 
to  enhance  the  range  resolution  to  a  level  that  far  exceeds 
the  limitation  discussed  earlier. 

Signal  processing  in  range-Doppler  imaging  can  be  con¬ 
sidered  as  a  two-dimensional  spectral  analysis  problem. 
Many  high-resolution  spectral  estimation  and  frequency 
estimation  methods,  i.e.,  modem  spectral  estimation 
methods,  have  been  studied.^  In  order  to  take  advantage 
of  the  superiority  of  these  methods,  an  accurate  signal 
model  and  a  sufficient  signal-to-noise  ratio  (SNR)  are 
required.  Furthermore,  modern  spectral  estimation 
methods  are  far  more  complex  compared  to  the  classical 
Fourier  method.  For  this  reason,  people  continue  to 
investigate  methods  that  combine  modern  high- 
resolution  techniques  with  the  classical  approach.  One  of 
the  hybrid  approaches^  is  relatively  simple  in  concept 
and  in  computation.  In  this  approach,  a  linear  prediction 
filtering  method  is  employed  to  extrapolate  the  data 
beyond  the  observation  window.  Then,  the  extrapolated 
data  undergoes  discrete  Fourier  transform  (DFT)  to 
attain  a  higher  frequency  resolution.  This  method  is  used 
to  improve  the  resolution  of  rotating  platform  imaging  to 
construct  a  two-dimensional  image  of  a  B-52  model  by 
processing  data  obtained  in  a  microwave  anechoic 
chamber.  Preliminary  results  indicate  that  this  superres¬ 
olution  range-Doppler  radar  imaging  method  is  less 
complex  in  computation  and  more  general  in  applica¬ 
bility.  In  the  area  of  modem  spectral  estimation  methods 
to  improve  radar  imaging  resolution  capability, 
Gabriel^’^  introduced  a  maximum  likelihood  filter 
method  and  a  linear  prediction  filter  method  to  achieve 
superresolution  range-Doppler  imaging.  These  methods 
involve  the  determination  of  the  inverse  of  matrices, 
which  requires  a  great  deal  of  computation.  Further¬ 
more,  references  4  and  5  only  provided  imaging  results 
of  simulated  data  and  did  not  verify  the  methods  with 
measured  data. 

In  this  paper,  the  use  of  linear  prediction  for  data 
extrapolation  to  achieve  superresolution  imaging  is 
described  first.  A  similar  method  is  then  used  to  obtain 
superresolution  of  cross-range  imaging.  Finally,  a  super¬ 
resolution  two-dimensional  image  of  a  B-52  model  is 
presented. 
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11.  Superresolntion  Range  Imaging 

The  high  range-resolution  capability  needed  in  range- 
Doppler  radar  imaging  can  be  realized  by  using  tech¬ 
niques  such  as  bandwidth  pulse  compression,  time- 
bandwidth  interchange,  or  frequency-stepping 
waveform.  The  method  described  here  is  also  applicable 
to  pulse  compression  and  time-bandwidth  interchange. 
The  difference  is  that  the  pre-treatment  is  slightly 
different.*’^ 

When  frequency-stepping  waveform  is  employed,  the 
radar  transmission  frequency  and  the  receiver  local- 
oscillation  frequency  vary  in  a  stepwise  manner.  Figure  1 
shows  a  target  at  a  distance  R  away  from  the  radar;  fl[r)  is 
the  range  function  of  the  reflectivity  of  the  target.  It  is 
also  known  as  the  range  image  of  the  target.  The  variable 
r  represents  the  distance  of  the  target  with  respect  to  the 
reference  point  o.  With  a  specific  viewing  angle,  the 
radar  transmits  a  series  of  K  pulses  at  a  frequency 

fk  =  fo  +  k8f,k  =  0,  1, . ,K-  1.  (3) 

The  return  from  the  target  after  compensating  for  the 


Radar  Target/^ 

Figure  1.  Target  Geometry 


phase  term  due  to  range  R  can  be  expressed  as  follows: 
5*=Aj/(r)  exp  {+/^/,r}4r, 

^=0,  1,  . A'-l.  (4) 


where  A  is  a  proportionality  factor  which  is  dependent 
upon  factors  such  as  power  and  signal  transmission.  It  is 
also  assumed  the  pulse  width  covers  the  entire  target 
range.  Integration  is  carried  out  along  r  over  the  entire 
target  range.  Equation  (4)  implies  that  is  the  inverse 
Fourier  transform  of  fl[r).  Hence,  after  knowing  the  K 
vector  of  S^,  it  is  possible  to  estimate  f(r)  by  computing 
the  DFT  of  where  k  =  0,  1, ....,  K  -  1.  Obviously,  the 
range  resolution  of  such  a  target  image  is 

p,  =  C/2K6f  (5) 

This  is  consistent  with  equation  (1)  because  here  the 
effective  bandwidth  of  the  transmitted  signal  is  B  =  K8f. 

Let  us  assume  that  the  scattering  from  an  extended  target 
comes  from  a  finite  number  of  scattering  centers.^  For  an 
individual  scattering  center,  if  its  range  is  ro,  the 
observed  signal  in  the  frequency  domain  is  a  complex 
sine  function  of 


A/(ro)exp|+;-^  rj^ 

Therefore,  the  observed  signal  of  the  target  is  comprised 
of  a  series  of  complex  sine  functions  in  the  frequency 
domain.  It  can  be  accurately  predicted  by  a  linear 
prediction  filter  of  an  order  equal  to  or  greater  than  the 
number  of  complex  sine  functions.^  Although  this  state¬ 
ment  is  valid  only  when  there  is  no  noise  in  the  observed 
signal,  in  reality  it  is  still  possible  to  employ  linear 
prediction  to  reliably  extrapolate  the  data  outside  the 
observation  window  when  the  SNR  is  relatively  low. 

Assuming  that  actual  data  S|„k  =  0, 1,....,K- 1, has  been 
obtained  at  a  given  viewing  angle,  a  single-step  linear 
prediction  filter  of  order  P  is  used  to  extrapolate  the 
data.  The  filter  coefficient  is  ap,  p  =  1, ....,  P.  This  filter 
produces  a  forward  extrapolation  sample  as  follows: 

■S';r=S5'*_,ap,  h>K  (6) 

P-1 

The  backward  extrapolation  process  is  similar.  The 
difference  is  that  the  filtering  direction  is  reversed. 
Moreover,  the  filter  coefficients  are  complex  conjugates 
of  those  of  the  forward  filter.  In  other  words,  the  extrap¬ 
olation  takes  place  according  to  the  following: 


^<1  (7) 

P-  I 

Thus,  the  data  is  actually  measured  if  1  <  k  <  K  and  the 
data  is  extrapolated  when  k  >  K  and  k  <  1 .  The  forward 
extrapolation  length  is  equal  to  that  of  backward  extrap¬ 
olation.  The  Burg  maximum  entropy  (ME)  technique  is 
used  to  calculate  coefficients  ap,  p  =  1,  ....,  P.  A  vector 
with  N  data  points  (N  >  K)  can  be  obtained  with  the 
method  described  above.  This  N-point  extrapolated 
vector  replaces  the  original  K-point  vector  in  the  DFT 
computation.  The  range  image  thus  derived  can  be 
expected  to  have  a  much  better  range  resolution. 

Figure  2  shows  the  overall  range  image  of  a  B-52  model. 
This  is  the  result  of  computer  processing  based  on 
microwave  anechoic  chamber  data.  A  frequency  step¬ 
ping  technique  was  used  in  the  experiment.  The  number 
of  frequency  steps  is  K  =  32.  The  transmitted  signal  has 
an  effective  bandwidth  of  B  =  2.15  GHz.  Table  1  shows 
the  raw  data  obtained  in  the  microwave  anechoic 
chamber  at  a  specific  viewing  angle.  Figures  2(a)  and  (b) 
are  the  range  images  obtained  with  the  original  32  data 
points  and  with  128  extrapolated  points  with  a  Ham¬ 
ming  window,  respectively.  It  is  apparent  that  data 
extrapolation  by  linear  prediction,  in  conjunction  with 
conventional  Fourier  transformation,  can  provide  super¬ 
resolution  range  images. 
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Figure  2.  Range  Image  of  B-52  Model  at  a  Specific  Viewing  Angle,  (a)  with  32  original  points,  (b)  with  128  extnpo* 

lated  points 


Table  1.  Raw  Data  at  a  Specific  Viewing  Angle 


-.080  99 

.11520 

-.12826 

.10653 

-.17118 

.07572 

-.19275 

.00786 

-.169  40 

-.08068 

-.13269 

-.18262 

-.06945 

-.26023 

-.01409 

-.27528 

.048  96 

-.21143 

.08769 

-.08913 

.09758 

.06676 

.  .08207 

.19537 

.026  67 

.26618 

-.04690 

.25384 

-.11979 

.18449 

-.17199 

.06819 

-.184  57 

-.05550 

-.16452 

-.14741 

-.10523 

-.19206 

-.02174 

-.19634 

.064  16 

-.17319 

.13968 

-.12860 

.19338 

-.08315 

.21725 

-.03815 

.214  41 

.01149 

.16985 

.06559 

.11442 

.13303 

.05429 

.17366 

-.015  34 

.17477 

-.07138 

.13274 

-.11886 

.03660 

-.14958 

-.06196 

III.  Snperresolution  Cross-Range  Imaging 

In  small-rotation-angle  imaging  with  the  FFT  range- 
Doppler  process,  range  and  cross  range  are  processed 
separately.  As  described  in  Section  II,  we  can  obtain  a 
range  image  at  every  viewing  angle.  At  M  equally  spaced 
viewing  angles,  M  range  images  of  the  target  can  be 
obtained.  Thus,  a  data  vector  of  M  points  can  be 
obtained  with  respect  to  each  range  unit.  The  cross-range 
image  of  a  given  range  unit,  i.e.,  the  reflectivity  of  the 
target  in  that  range  unit  as  a  function  of  the  cross  range, 
can  be  derived  from  the  DFT  of  this  M-point  vector. 
When  it  is  applied  to  a  sample  of  M  actually  measured 
data  points,  the  cross-range  resolution  is 

-  X/2M8e  (8) 

where  86  is  the  spacing  of  two  neighboring  viewing 
angles.  Since  A0  -  M80,  equation  (8)  is  consistent  with 
equation  (2). 


On  the  basis  of  the  same  argument  to  attain  super  raiife 
resolution  discussed  in  Section  II,  superresolution  crosh- 
range  images  can  be  obtained  by  using  linear  prediction 
for  data  extrapolation  in  conjunction  with  conventional 
Fourier  transformation.  Figure  3  shows  the  overall  cross- 
range  image  of  the  B-52  model  at  a  range  unit  derived 
from  measured  data.  The  total  angle  of  rotation  of  the 
airplane  model  in  the  imaging  process  is  AO  -  1 S*  and  the 
radar  wavelength  is  A.  ••  2  cm.  Table  2  lists  the  32  original 
data  points  corresponding  to  the  range  unit.  Figurm  3(a) 
and  (b)  show  the  result  obtained  from  the  original  32 
data  points  and  that  derived  from  the  128  extrapolated 
points  with  a  Hamming  window,  respectively.  It  is 
apparent  that  using  linear  prediction  for  data  extrapola¬ 
tion  can  also  significantly  improve  the  cross-range  reso¬ 
lution  of  range-Doppler  radar  imaging. 


:,JP4^STr»2r023 
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Figure  3.  Cross-Range  Image  of  B-52  Model,  (a)  with  32  original  points,  (b)  with  128  extrapolated  points 


Table  2.  Raw  Data  at  a  Given  Range  Unit 


t.391  *29  '  ’ 

.43871  ^ 

-.07534 

-.23074 

.50295 

-.80743 

.53297 

1  -.65919 

-.569  85 

-.42205 

-1.20443 

-.69226 

-.40890 

-.36944 

.21238 

.43367 

.303  96 

.82409 

.19226 

.46031 

.17158 

-.13421 

.11666 

-.61985 

-.697  28  -  - 

-.70557 

-.27174 

-.52000 

-.16635 

-.05365 

.10964 

.395  75 

.54619 

.62743 

.19806 

.82085 

-.28259 

.75538 

-.51553 

.319  49 

-.41673 

-.17068 

-.27398 

-.13106 

-.10458 

.21626 

.16420 

.479  72 

.42178 

.43015 

.29670 

.23022 

.  .-..1 

-.15925 

.01776 

-.47091 

-is?  82  .  -  . 

-.43817 

-.51430 

-.18886 

-.53399 

-.04538 

-.26843 

-.02514 

^lYh  Sui^rresqlution  Two-Dimensional  Imaging 

The  range  resolution  and  cross-range  resolution  in  two- 
dimensional  innge-Doppler  ima^ng  can  both  be  improved 
%fth  the  ttiethod  described  above.  The  results  are  shown  in 
Figuro  4,  Using  the  data  shown  in  two  previous  sections,  a 
, two-dimensional  radar  image  of  the  aiiplane  model  was 
^pbtained  by  processing  the  range  dimension  first  for  every 
/viemng  angle  arid  then  processing  the  cross  range  for  every 
'  tarige  uriii.  Figures  4(a)  and  (b)  show  the  images  of  the  B-52 
model  processed  based  on  the  original  32  x  32  data  matrix 
^^nd  the  extrapolated  128  x  128  matrix,  respectively.  A 
.  ^9pmparison  shows  that  our  approach  significantly  improved 
;  range  resolution  and  cross-range  resolution  of  the  image. 
(Cofrtprir^d  to  the  photograph  of  the  B-52  model  shown  in 
Figure  5,  the  two  engines  Under  the  wing  were  clearly  visible 
in.  the  superresolution  two-dimensional  image  shown  in 


Figure  4(b). 


Figure  4.  Two-Dimensional  Image  of  B-52  Model,  (a) 
32  X  32  true  data  matrix,  (b)  128  x  128  extrapolated 
data  matrix 


Figure  5.  Pbotograpb  of  B-52  Model 
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The  microwave  anechoic  chamber  data  used  has  a  very  high 
SNR.  It  almost  can  be  treated  as  a  noise-free  case.  In  order 
to  investigate  the  effect  of  SNR  on  the  superresolution 
imaging  technique  described,  white  Gaussian  noise  was 
added  to  the  microwave  anechoic  chamber  data.  A  great 
deal  of  results  indicate  that  this  method  can  still  produce 
better  superresolution  two-dimensional  images  even  with  a 
low  SNR.  Figure  6  shows  the  results  with  an  SNR  of  0  dB. 


Figure  6.  Two-Dimensional  Images  of  B-52  Model  With 
SNR  =  0  dB.  (a)  32  x  32  true  data  matrix,  (b)  128  x 
128  extrapolated  data  matrix 


V.  Conclusions 

Our  preliminary  results  indicate  that  superresolution 
range-Doppler  radar  imaging  using  data  extrapolation 
via  linear  prediction  is  a  relatively  simple  and  practical 
superresolution  imaging  method.  It  is  capable  of  lifting 
the  resolution  limitation  of  conventional  Fourier  image 
processing.  This  method  can  produce  higher  resolution 
images  with  identical  effective  transmission  bandwidth 
and  total  angle  of  rotation  during  imaging.  Or,  it  can 
produce  an  image  of  similar  resolution  with  a  smaller 
bandwidth  and  total  rotation  angle. 

During  the  process  of  revising  this  paper,  the  authors 
found  that  Li*  et  al.  used  a  similar  method  to  improve 
the  resolution  in  one-dimensional  range  imaging  by 
extrapolating  data  in  different  isolated  segments.  This 
work  has  presented  a  method  for  two-dimensional  super- 
resolution  range-Doppler  imaging.  In  addition,  the  effect 
of  SNR  on  this  method  has  also  been  studied.  Even  with 
a  low  SNR,  this  method  was  still  found  to  be  capable  of 
producing  better  superresolution  two-dimensional 
images. 

The  authors  wish  to  express  their  gratitude  to  Professors 
Liu  Yongtan  [0491  3057  0962],  Cao  Zhidao  [2580  1807 
6670],  and  Meng  Xiande  [1322  2009  1795]  of  Harbin 
Institute  of  Technology  and  to  colleagues  at  the  207th 
Institute  of  the  Ministry  of  Aerospace  Industry  for 
providing  the  microwave  anechoic  chamber  data. 
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[Text]  Abstract 

The  imaging  of  targets  rotating  in  a  small  azimuth  angle 
range  is  discussed.  A  unified  image  reconstruction 
scheme  suitable  for  bistatic,  monostatic,  near-field  and 
far-field  imaging  geometries  is  derived.  Super-resolution 
imaging  techniques  are  investigated.  A  new  imaging 
reconstruction  algorithm,  combining  spectral  extrapola¬ 
tion  and  coarse  “shape”  information,  is  presented. 
Finally,  experimental  results  in  microwave  anechoic 
chambers  are  described. 

Key  Words:  microwave  imaging,  radar  targets,  imaging 
processing. 

I.  Introduction 

Microwave  imaging  of  rotating  targets  has  been  widely 
used  in  many  areas,  such  as:  (1)  Target  diagnosis  and 
identification.  Rotating  target  imaging  is  a  basic  tool  for 
studying  the  electromagnetic  scattering  of  targets.  It 
plays  an  important  role  in  the  study  of  electromagnetic 
scattering  mechanism  and  RCS  reduction  for  stealth 
research  and  target  identification.  (2)  SAR  and  ISAR 
signal  processing.  SAR  or  ISAR  imaging  with  motion 
compensation  is  equivalent  to  rotating  target  imaging. 
(3)  Non-destructive  testing.  The  principle  of  rotating 
target  imaging  is  also  widely  used  in  non-invasive  and 
non-destructive  testing.  Therefore,  in-depth  study  of 
rotating  target  imaging  and  related  processing  tech¬ 
niques  is  a  basic  research  project  with  a  wide  range  of 
applications. 

D.  L.  Mensa^  once  pointed  out  that  imaging  of  targets 
rotating  over  360°  could  achieve  extremely  high  nominal 
two-dimensional  resolution  (a  fraction  of  a  wavelength 
of  the  radar).  The  condition  to  attain  such  a  high 


resolution  is  that  all  the  scatterers  on  the  target  must  be 
isotropic.  Furthermore,  the  scatterers  must  not  move  in 
order  to  obtain  a  uniform  synthetic  radar  aperture  over 
360°.  However,  the  following  two  reasons  make  rotation 
imaging  in  small  azimuth  angle  range  more  practical.  (1) 
Real  radar  targets  often  cannot  satisfy  the  requirement 
that  the  scatterers  be  isotropic  and  stationary.  Hence,  it 
is  not  possible  to  attain  the  high  theoretical  nominal 
resolution.  Furthermore,  the  amount  of  work  involved  in 
data  acquisition  and  processing  for  360°  rotating  target 
imaging  is  tremendous.  (2)  As  radar  signal  technology 
advances,  high  radial  resolution  of  the  target  can  be 
obtained  with  a  wideband  signal  alone.  It  usually  takes  a 
small  azimuth  angle  to  achieve  high  transverse  resolu¬ 
tion.  Therefore,  imaging  of  targets  in  a  small  azimuth 
angle  range  has  been  developed  for  a  variety  of  applica¬ 
tions  in  the  past  two  decades. 


II,  Image  Reconstruction  Algorithm 

The  bistatic  geometry  for  rigid  target  imaging  is  shown  in 
Figure  1.  The  target  coordinate  is  (x-y)  and  the  radar 
coordinate  is  (u-v). 


X  =  u  cos  0  -  V  sin  8, 
0 


y  =  u  sin  0  +  V  cos 
(1) 


or 

u  =  X  cos  0  +  y  sin  0,  v  =  -  x  sin  0  +  y  cos 

0  (2) 


antenna  antenna 

Figure  1.  Geometry  for  Microwave  Imaging 
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The  two-way  range  between  a  point  on  the  target  (x,y) 
and  the  transmitter  and  receiver  is 

R.=L(.Ro+vy+(.d+uyy^ 

=  {Rl  +y^+d^+2r[IioSm  (fi-O) 

+dcos(f>-0)jy"  (3) 

Rr=[(Ra+vy+(d-uyy'^ 

=  {Ri  +  y+d^+2r[JiaSin  (P-0) 

-dcos(p-0)2}*'2  (4) 

and 

R=R,+Rr 

=  {RS  +  y-td^+2r[Jiism  (P-0) 

■ficos  (P  — 0)]}’'^ 

+  {^o  +  rHi*-f2rC/2osin  (P—0) 

.  -icoS(P-0)]}''* 

From  the  above  formulas,  after  frequency  mixing,  cali¬ 
bration  and  filtering,  the  return  can  be  derived  as  fol¬ 
lows: 

=jj  J(r, <p) exp {-j:iK(.R,-\-Rr 

-2Ra'>}Tdrdf  (6) 

From  equation  (6),  the  return  is  the  correlation  integral 
of  the  scattering  signals  along  various  equidistance 
lines;  fl[r,(p)  is  the  scattering  function  of  the  target.  In 
two-dimensional  imaging  of  rotating  targets,  the  spec¬ 
tral  domain  signal  of  fl(r,(p)  is  measured  in  the  ring- 
shaped  spectral  domain  shown  in  Figure  2.  The  inte¬ 
gration  for  image  reconstruction  is  also  performed  in 
this  domain.  The  image  reconstruction  equation 
expressed  in  polar  coordinate  is  as  follows: 

/(r,|>)  Xexp{;--rJC(/e , 

--2R;>}dKde 


Figure  2.  Polar  Coordinate  Data  Format  in  Imaging  of 
Rotating  Target 


Various  reconstruction  methods  can  be  designed  based 
on  equation  (7),  including  the  two-dimensional  Fast 
Fourier  Transform  (2-D  FFT),  filtered  back-projection 
(FBP)  and  Fourier  domain  extrapolation. 

Experimental  and  theoretical  studies  indicate  that  the 
FBP  method  not  only  is  more  flexible  but  also  produces 
higher  quality  reconstructed  images  in  two-dimensional 
imaging  of  rotating  targets.  It  is  suitable  for  anechoic 
chamber  imaging.  This  method  is  described  as  follows: 

In  equation  (7),  if  the  integration  of  K  is  treated  as  the 
inverse  Fourier  transform  of  variable  ^  - 

Ro)/2,  then  the  reconstructed  ftcaret](r,(p)  becomes 

P,(.l,)=l‘^^~’^'\K  +  Ki\F(K+K„e) 
exp{}2‘-rKl,}dK 

fir, 9)  ,il  .'>txp{j2^Kil  ,}ie 

l,  =  (.R,  +  Rr-Ro')/2  (8) 


In  actual  computation,  the  first  formula  (radial  focusing) 
is  computed  by  FFT  and  then  the  second  formula  is  used 
to  compute  back-projection  (transverse  focusing).  It  is 
very  similar  to  the  FBP  technique  in  tomographic 
imaging.  Hence,  it  is  also  called  filtered  back- 
projection.^’^ 

The  third  formula  in  equation  (8)  describes  the  path  of 
back-projection.  From  analytical  geometry,  the  path  of 
back-projection  described  by  that  formula  follows  an 
elliptical  trajectory.  It  is  identical  to  the  integration  path 
of  radar  projection  measurement,  i.e., 

ir  =  {:i?S  +  H+d*-|-2r(J?osin  (P-0) 

+  dCOS 

+  11/22  + H-|-i*-|-2r(i?„ sin  (P-0) 

- d  cos  (P  -  0)  )  2Ro}  /2  (  9  ) 

If  it  is  a  monostatic  measurement,  then  in  Figure  1  d  ==  0 
and  equation  (9)  becomes 

le  -  [R2  +  r^  +  2r  Ro  sin  ((p  -  e)]''^  ^  (jO) 

This  describes  a  family  of  circles.  In  other  words,  in 
monostatic  measurements,  the  projection  integration 
path  is  a  circular  arc  and  the  back-projection  path  in 
image  reconstruction  is  also  circular. 
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Furthermore,  when  Rq  >>  r,  i.e.,  the  far-field  condition 
is  met,  then 

Ic  «  r  sin  ((p  -  6)  (11) 

Equation  (11)  describes  a  straight  line.  It  is  suitable  in 
far-field  imaging  (plane  wave  front). 

The  derivation  of  equation  (8)  does  not  require  the  radar 
to  emit  plane  electromagnetic  wave  fronts.  Therefore, 
the  resultant  reconstruction  method  is  applicable  to  a 
variety  of  imaging  conditions  such  as  monostatic, 
bistatic,  far-field  and  near-field.  Incidentally,  because 
the  target  spectrum  obtained  in  microwave  imaging  of  a 
rotating  target  is  located  in  a  ring-shaped  domain,  the 
expansion  function  of  the  imaging  system  has  very  high 
side  lobes.  Therefore,  usually  a  smoothing  window  (at 
small  azimuth  angle)  or  variable  trajectory  filter  (360° 
imaging)^  is  required  in  the  actual  reconstruction  com¬ 
putation  to  suppress  the  side-lobe  voltage  in  the  recon¬ 
structed  image  to  improve  image  quality. 


III.  Super-Resolution  Imaging  Technique 

It  is  well  known  that  linear  prediction  is  widely  used  in 
the  estimation  of  high-resolution  power  spectrum  of  a 
discrete  time  sequence,  especially  when  the  sequence  is 
very  short.  By  means  of  extrapolating  the  sequence,  it  is 
possible  to  achieve  higher  resolution  than  that  of  a 
conventional  Fourier  spectrum.  Since  in  microwave 
imaging  the  two-dimensional  resolution  is  determined 
by  the  frequency  spectrum  window  and  azimuth 
window,  when  the  radar  bandwidth  is  narrow  or  the 
azimuth  angle  is  small,  the  resolution  deteriorates  rap¬ 
idly.  In  order  to  obtain  a  high-resolution  microwave 
image  with  a  small  amount  of  data,  linear  prediction  is 
used  to  investigate  the  super-resolution  imaging  of 
microwave  data. 

Based  on  the  Burg  algorithm,^*^  a  sequence  with  M 
specimens  [xj,  X2,  X3,  ...,  x^]  is  extrapolated  by  linear 
prediction  to  have  N  (N  >  M)  specimens.  The  predic¬ 
tion  parameters  can  be  obtained  by  minimizing  the  sum 
of  errors  in  energy,  Ep,  of  forward  and  backward  predic¬ 
tion,  i.e.. 


(12) 


where  Epf  is  the  forward  prediction  error  of  the  pth  order 
model: 


Epf 


p 


♦  -0 


(13) 


and  Epb  is  the  backward  prediction  error  of  the  pth  order 
model: 

E  ph  it  (14) 


Here,  apj^  is  the  linear  prediction  parameter  and  ♦  repre¬ 
sents  a  complex  conjugate  operation. 

The  linear  prediction  equation  for  extrapolation  with  M 
specimens  is  as  follows: 

A  ^ 

/>0  (15) 

^-1 

*-/=S®**^-/**  ;>o,  (16) 


where  [caret]  represents  the  estimated  value. 


Figure  3  shows  schematic  diagrams  of  one-  and  two- 
dimensional  extrapolation  by  linear  prediction  with  par¬ 
tial  data  in  the  target  spectrum.  After  spectral  expansion, 
the  image  can  be  reconstructed  with  the  resultant  data 
using  the  method  described  in  the  last  section  to  obtain 
a  super-resolution  image  of  the  target. 


Azimuth 

spectral 

expansion 


Freauency 

spectral 

expansion 


(a)  One*«dimensional  spectral 
expansion 


VSH  Measured 
□  Predicted 


(b)  Two-dimensional  spectral 
expansion 

Figure  3.  image  Spectrum  Expansion  by  Linear 
Prediction 


In  image  reconstruction  with  extrapolated  data  by  linear 
prediction,  it  was  found  that  although  the  resolution  is 
enhanced,  the  side-lobe  voltage  of  the  image  is  also 
increased  due  to  prediction  errors.  In  order  to  suppress 
image  side  lobes,  an  algorithm  is  developed  to  combine 
linear  prediction  with  the  rough  “shape”  of  the  target. 
The  “shape”  is  derived  from  the  original  data.  Specifi¬ 
cally,  both  original  data  and  extrapolated  data  are  pro¬ 
cessed  in  image  reconstruction.  In  radial  focusing  with 
FBP,  two  sets  of  range  distribution  plots  are  obtained. 
The  range  distribution  plot  obtained  from  the  extrapo¬ 
lated  data  is  compared  to  that  derived  from  the  original 
data  to  eliminate  the  side  lobes.  After  both  images  are 
reconstructed,  comparison  can  be  made  again  to  further 
eliminate  any  side  lobes  other  than  the  “shape”  of  the 
target  in  the  original  image.  This  new  algorithm  can 
drastically  lower  the  side-lobe  voltage  of  the  image. 
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Figure  4  shows  the  flowchart  of  this  method.  The  results 
will  be  presented  in  the  next  section. 


i^^reauencv  speetriiin  expanslonj 
iAy.-iTOuty  angle  spec^ryffl  gxpanalflj 


ilmage  recoT^structionj 
target  i^aje] 

Figure  4.  Flowchart  of  the  New  Image  Reconstruction 
Algorithm 


lY.  Anechoic  Chamber  Experiments  and  Results 

The  entire  experiment  was  done  on  the  microwave 
measurement  and  imaging  system  of  an  HP8510  net¬ 
work  analyzer.  The  major  specifications  of  the  system 
are  shown  in  Table  1. 


Table  !•  Major  Specifications  of  the  Microwave  Imaging 
System 


2-4,  4-8,  8-18  GHz 
1  X  lO'^/day 


Sweep  range 
Frequency  stability 
Output  power 
Antenna  gain 

Rotating  platform  accuracy 

Minimum  step  of  rotating  plat¬ 
form 

Amplitude  measurement  accu¬ 
racy 

Sweep  linearity 
Sweep  mode 
Receiver  sensitivity 

Dynamic  range 

Phase  measurement  accuracy 


20  dBm 
23  dB 
0.05' 

o.r 

-I-/-  0.5  dB 

1%  in  sweeping  mode 

Stepping,  sweeping 

-85  dBm  (unaveraged) 

-103  dBm  (averaged  over  64 
times) 

75  dB  (unaveraged) 

93  dB  (averaged  over  64  times) 
+/-  3' 


During  the  imaging  process,  a  microcomputer  controls 
the  transmitter,  receiver  and  azimuth  angle.  The  ampli¬ 
tudes  and  phases  of  the  wide-band  scattered  returns  from 
the  target  are  collected  at  different  azimuth  an^es.  The 
data  gathered  is  also  pre-processed  by  the  microcom¬ 
puter. 


The  pre-processing,  including  software  range  gate  tech¬ 
nology  for  background  cancellation,  data  scaling  and 
clutter  filtering  outside  the  target  range,  can  be  expressed 
by  the  following  equation: 


n — p 


(17) 


where  E  is  the  accurate  amplitude  phase  of  the  cali¬ 
brator,  T  is  the  frequency-domain  measurement  of  the 
target,  Tj  is  the  intrinsic  scattering  of  the  target,  S  is  the 
frequency-domain  measurement  of  the  calibrator,  Sj  is 
the  inherent  scattering  of  the  calibrator,  is  the  back¬ 
ground  during  target  measurement,  is  the  background 
during  scalar  measurement,  and  script  F[  ]  represents  the 
filtering  process. 

The  RCS  amplitude  of  the  pre-processed  imaging  data 
can  be  expressed  in  dB  referenced  to  a  square  meter 
(dBsm).  The  phase  is  centered  around  the  rotating  plat¬ 
form.  The  data  is  then  used  for  image  reconstruction  and 
image  output  by  a  microcomputer. 

Figure  5  shows  the  resultant  image  of  a  reduced  B52 
model  in  an  anechoic  chamber.  The  experimental  con¬ 
ditions  are:  center  frequency  11.5  GHz,  radar  bandwidth 
7.0  GHz,  frequency  step  35  MHz,  target  rotation  angle 
20®,  angular  sampling  interval  0.5®,  antenna  polarization 
state  H-H,  test  range  12.5  m.  A  Hamming  window  was 
added  in  image  reconstruction.  The  specifications 
attained  are:  radial  resolution  approximately  4  cm, 
transverse  resolution  7  cm,  and  image  side-lobe  voltage 
less  than  -40  dB.  The  dynamic  range  of  the  imaging 
system  is  actually  limited  by  the  side-lobe  voltage.  The 
dotted  line  in  the  figure  represents  the  actual  silhouette 
of  the  reduced  model.  Its  scattering  center  and  the 
contour  lines  obtained  from  imaging  are  in  excellent 
agreement. 

Figures  6  and  7  show  the  reconstructed  images  derived 
from  the  raw  data  and  the  two-dimensional  extrapolated 
data,  respectively.  Some  data  were  acquired  with  a 
reduced  bandwidth  of  3.5  GHz  and  smaller  azimuth 
angle  of  15®.  Comparing  Figure  7  to  Figure  5,  we  find 
that  the  resolution  of  the  reconstructed  image  from  the 
raw  data  and  that  of  the  reconstructed  image  based  on 
extrapolated  data  are  essentially  the  same.  However, 
there  is  some  error  in  the  amplitude  of  the  scattering 
center  (approximately  1-2  dB). 


V.  Conclusions 

(1)  Using  a  7  GHz  wide-band  frequency  spectrum  and 
20®  rotation  angle  aperture  imaging,  it  is  possible  to 
obtain  a  high-resolution  microwave  image  of  the  order  of 
a  centimeter. 

(2)  The  FBP  algorithm  derived  is  applicable  to  a  variety 
of  imaging  geometries  and  conditions  including  bistatic, 
monostatic,  near-field  and  far-field. 

(3)  By  combining  spectral  expansion  and  target  rough- 
shape  data,  super-resolution  imaging  can  be  attained. 
Consequently,  the  amount  of  imaging  data  required  is 
drastically  reduced.  At  3.5  GHz  and  15®  azimuth 
window,  the  resolution  is  almost  as  high  as  at  7  GHz  and 
20®  azimuth  window,  and  side-lobe  voltage  is  better  than 
-40  dB.  However,  there  is  some  error  in  amplitude. 


dBnn 


Transverse  distance/cm 


Transverse  distance/cm 

(a)  Three-dimensional  diagram  (b)  Contour  diagram 

Figure  5.  Imaging  of  a  B52  Model 


Transverse  distance/cm 

(a)  Three-dimensional  diagram 
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Figure  6.  Image  Reconstructed  From  Raw  Data 
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Transverse  distance/cm 
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Figure  7.  Image  Reconstructed  From  Spectral  Expanded  Data 
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Finally,  in  addition  to  spectral  expansion,  the  quality  of 
the  image  may  be  further  improved  by  using  polarization 
information  and  target  symmetry.  This  effort  is  in 
progress. 
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[Excerpts]  Abstract 

In  this  work,  the  scattering  of  a  large  planar  slot-array 
antenna  is  studied  using  the  momentum  method,  phys¬ 
ical  optics  method  and  singularity  elimination  tech¬ 
nique.  The  radar  cross  section  (RCS)  of  a  typical  array 
was  calculated  and  verified  experimentally.  It  is  found 
that  the  calculated  result  is  in  good  agreement  with  the 
experimental  data.  The  computations  also  show  that  the 
contribution  of  the  plane  to  the  scattering  field  (structure 
term)  is  far  greater  than  that  of  the  slot  (mode  term). 

Key  Words:  planar  slot-array  antenna,  radar  cross  sec¬ 
tion,  scattering. 

I.  Introduction 

The  large  planar  slot-array  antenna  is  one  of  the  com¬ 
monly  used  airborne  antennas.  However,  it  produces  a 
very  strong  RCS  contribution  in  the  direction  of  the  nose 
of  the  aircraft.  In  order  to  improve  the  stealth  character¬ 
istics  of  the  aircraft,  it  is  necessary  to  understand  the 
scattering  mechanism  of  this  type  of  antenna  and  ana¬ 
lyze  and  calculate  its  RCS  contribution.  To  date,  most  of 
the  studies  on  the  coupling  of  incident  wave  and  slot- 
array  antenna,  or  scattering  computation,  involved  the 
combination  of  a  single  slot  and  a  cavity. There  is  no 
information  available  on  the  scattering  mechanism  of  a 
large  planar  slot-array  antenna  and  the  detailed  analysis 
of  its  RCS.  This  work  provides  a  theoretical  analysis  of 
this  type  of  antenna  by  rigorously  taking  the  interaction 
between  elements  into  consideration  using  momentum 
method,  physical  optics  method  and  singularity  elimina¬ 
tion  technique.  Moreover,  the  computing  time  is  drasti¬ 
cally  reduced  through  the  use  of  symmetry  technique, 
Toeplitz  property  and  an  excellent  inverse  technique. 

IL  Antenna  Array  Analysis 

Figures  1  and  2  show  the  antenna  array.  Let  us  assume 
that  the  incident  direction  is  on  the  yoz  plane.  The 
number  of  slot  vibrators  along  the  x-axis  over  a  length  of 
21  is  M,  the  distance  between  slots  is  Ax,  the  slot  width  is 
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W,  the  number  of  vibrators  long  the  y-axis  is  N  and  the 
spacing  is  Ay.  The  contribution  of  the  entire  antenna 
array  to  the  RCS  is  the  phase  sum  of  two  terms,  i.e.,  the 
mode  term  and  structure  term.^ 


Figure  1.  Array  Antenna 


Zone  a :  Zone  b 

l^onductor 

Incident 
direction  i 

! 

yl 

Figure  2,  Side  View  of  Array 
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Figure  3.  Equivalent  Zones 
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the  axial  direction.  Hence,  the  resultant  far  scattering 
field  essentially  varies  according  to  the  pattern  of  the 
plate  structure  term  in  the  range  of  0-60“.  Therefore,  it  is 
quantitatively  demonstrated  that  structural  scattering 
from  the  plate  is  the  dominating  term  in  the  mechanism 
of  a  slot-array  antenna.  In  order  to  effectively  reduce  the 
RCS  of  a  slot-array  antenna,  it  is  necessary  to  minimize 
the  contribution  of  the  plate.  For  example,  a  traveling 
wave  slot-array  antenna  (non-resonant)  may  be  used  to 
offset  the  peak  of  its  emission  from  its  peak  direction  of 
scattering  to  lower  the  RCS  at  the  nose  of  the  aircraft.  In 
addition,  a  conformal  curved  surface,  such  as  a  cylin¬ 
drical,  conical,  spherical,  ellipsoidal  or  paraboloidal  sur¬ 
face,  can  be  used  instead  of  a  flat  plate  to  significantly 
reduce  the  contribution  of  the  structure  term  to  the 
RCS.® 


1.  Computation  of  Mode  Term 

The  plate  can  be  considered  infinite.  Its  conductive  wall 
divides  the  space  into  zones  a  and  b.  Using  the  equiva¬ 
lence  theorem  [See  Reference  5,  Sec.  3-5],  the  problem 
can  be  split  into  two,  as  shown  in  Figure  3.  When  the 
aperture  is  blocked  by  the  conductive  plate  and  replaced 
with  an  equivalent  magnetic  flux  M,  [I  =  (m-l)N  +  n,  m 
=  1,2, ...,  M,  n  =  1,  2, ...,  N],  the  field  in  zone  a  remains 
unchanged.  When  the  aperture  is  blocked  by  the  conduc¬ 
tive  plate  and  replaced  by  an  equivalent  magnetic  flux 
-Mj,  the  field  in  zone  b  remains  constant  as  well,  [passage 
omitted] 

III.  Results  and  Analysis 

As  an  example,  computation  was  done  on  a  4  x  8 
half-wavelength  slot-array  antenna  with  X  «  3.2  cm.  The 
slot  width  is  W  -  0.0 125^  The  end-point  spacing  in  the 
X  direction  is  a  half-wavelength  and  the  slot  spacing  in 
the  y  direction  is  also  a  half-wavelength.  The  plate  size  is 
4.0625^  X  lOX.  Figure  4  shows  the  RCS  of  the  antenna  as 
a  function  of  observation  angle  0.  The  unit  of  RCS  in  the 
figure  is  dBX^,  From  the  figure  one  can  see  that  the 
mode-scattering  term  generated  by  the  excitation  of  the 
incident  wave  at  the  slot  is  far  less  than  the  structure 
term  when  the  mutual  coupling  between  slots  is  taken 
into  consideration.  For  instance,  it  is  21.340  dB  less  at 


Figure  4.  RCS  of  Planar  Slot-Array  Antenna  vs.  Obser¬ 
vation  Angle 


An  experimental  antenna  was  made  based  on  the  param¬ 
eters  of  the  model  used  in  the  computation.  The  RCS 
was  measured  with  an  outdoor  quasi-monostatic  RCS 
measurement  system  (see  Figure  5).  The  measured  result 
is  compared  to  the  theoretical  value,  and  they  are  found 
to  be  in  good  agreement.  The  rotation  an^e  a  corre¬ 
sponding  to  the  maximum  RCS  value  in  Figure  5  should 
correspond  to  0  *  0“  in  Figure  4. 
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Platform  rotation 
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Figure  5.  Measured  RCS  of  Planar  Slot-Array  Antenna 


lY.  Conclusions 

Theoretical  analysis  and  numerical  computation  of  the 
RCS  of  a  planar  slot-array  antenna  was  carried  out  using 
the  momentum  method,  physical  optics  method  and 
singularity  elimination  technique.  It  was  found  that  the 
structure  term  plays  a  dominant  role  in  the  scattering 
mechanism.  Viable  technical  approaches  of  engineering 
value  are  presented  to  reduce  RCS. 
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[Excerpts]  Abstract 

The  backscattering  of  an  electromagnetic  missile  with  its 
energy  decaying  as  r"®  by  an  ideal  metallic  elliptical 
cylinder  and  a  sphere  is  described.  The  study  shows  that 
the  rates  of  energy  decay  for  the  backscattered  fields 
from  the  elliptical  cylinder  and  sphere  are  proportional 
to  and  respectively. 

Key  Words:  electromagnetic  missile,  backscattering, 
backscattered  field. 

I.  Introduction 

The  concept  of  an  electromagnetic  missile  was  first 
presented  by  T.  T.  Wu  in  1985.  It  is  a  transient  electro¬ 
magnetic  wave  of  finite  energy  that  does  not  decay 
according  to  as  a  conventional  spherical  wave  in  the 
far  field.  Instead,  it  decays  slowly  at  an  arbitrary  rate. 
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This  excellent  slow  decay  characteristic  of  the  electromag¬ 
netic  missile  can  be  used  to  increase  the  effective  range  and 
improve  the  performance  of  radars.*  Properties  such  as 
short  duration  (picoseconds)  and  high  peak  power  (giga¬ 
watts)  can  also  be  used  to  enhance  the  detection  capability 
of  the  radar  for  low-RCS  or  stealth  targets.  Therefore,  there 
is  a  need  to  have  a  profound  understanding  of  the  genera¬ 
tion  and  propagation  of  electromagnetic  missiles  from  both 
theoretical  and  practical  viewpoints. 

In  the  past  few  years,  a  great  deal  of  theoretical  and 
experimental  work  has  been  done  on  electromagnetic 
missiles.  Several  electromagnetic  missile  launchers  were 
constructed.  This  demonstrated  that  it  is  not  that  diffi¬ 
cult  to  make  electromagnetic  missiles  a  reality.  Besides 
constructing  a  variety  of  electromagnetic  missile 
launchers,  there  is  another  research  topic  of  considerable 
significance,  i.e.,  the  study  of  scattering  characteristics  of 
electromagnetic  missiles  from  different  target  shapes. 

It  is  well  known  that  the  energy  of  a  conventional  electro¬ 
magnetic  wave  decays  according  to  r"*  in  two-dimensional 
space.  Hence,  the  energy  of  scattered  field  decays  as  a 
function  of  r^.  In  three-dimensional  space,  the  energy  of  a 
conventional  electromagnetic  field  decays  as  a  function  of 
and  that  of  the  corresponding  scattered  field  decays 
according  to  r^.  Recently,  Myers  and  T.  T.  Wu  have 
studied  the  backscattering  of  electromagnetic  missiles 
from  a  plate.  The  results  show  that  the  energy  of  the 
backscattered  field  of  an  electromagnetic  missile  decaying 
as  a  function  of  r®  from  a  flat  plate  still  decays  according 
to  r®  (0  <  e  <  1).  This  is  much  slower  than  the  decay  of 
the  energy  of  a  conventional  scattered  field. 

This  study  involves  the  backscattering  of  electromag¬ 
netic  missiles  from  ideal  elliptical  cylindrical  and  spher¬ 
ical  conductors  of  infinite  length.  The  results  indicate 
that  the  backscattered  fields  of  the  elliptical  cylinder  and 
sphere  decay  according  to  r"^*"^®^  jK2+e)^  respectively. 
They  are  considerably  faster  than  the  decay  of  the 
backscattered  field  from  a  flat  plate.  However,  they  are 
slower  than  the  decays  of  conventional  two-  and  three- 
dimensional  scattered  fields,  which  are  proportional  to 
and  r*^,  respectively,  [passage  omitted] 


IV.  Conclusions 

This  paper  describes  the  backscattering  of  an  electro¬ 
magnetic  missile  decaying  as  a  function  of  r®  from  an 
elliptical  cylinder  conductor  and  a  sphere  conductor.  An 
analysis  shows  that  their  scattered  fields  are  proportional 
to  and  respectively.  This  conclusion  is  an 
important  reference  point  in  counter-stealth  technology. 
As  far  as  counter-stealth  is  concerned,  an  electromag¬ 
netic  missile  is  an  excellent  electromagnetic  pulse.  When 
used  in  impulse  radars,  it  can  effectively  counter  the 
effects  of  [radar]  absorbent  materials  (RAM)  and  profile 
stealth  technology.  The  signal-to-noise  ratio  of  the  radar 
can  be  improved  and  the  effectiveness  of  the  radar  can 


be  enhanced.  It  is  of  great  significance  to  further  inves¬ 
tigate  the  theop^  and  practice  of  backscattering  of  elec¬ 
tromagnetic  missiles  from  other  targets.  There  is  a  great 
deal  of  work  still  to  be  done  in  this  field. 
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